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Abstract. We present new parameterizations of vector and axial nucleon form factors. We 
maintain an excellent descriptions of the form factors at low momentum transfers (Q^), where 
the spatial structure of the nucleon is important, and use the Nachtman scaling variable ^ to 
relate elastic and inelastic form factors and impose quark-hadron duality constraints at high 
where the quark structure dominates. We use the new vector form factors to re-extract 
updated values of the axial form factor from j/^ experiments on deuterium. We obtain an 
updated world average value from ff^d, V^Y{ and pion electroproduction experiments of Ma 
= 1.014 ± 0.014 GeV/c? . Our parameterizations are useful in modeling v interactions at low 
energies (e.g. for oscillations experiments). The predictions for high can be tested in 
the next generation electron and Ufj, scattering experiments. (Presented by A. Bodek at the 
European Physical Society Meeting, EPS2007, Manchester, England, July 2007). 



The nucleoli vector and axial elastic form factors have been measured for more than 50 years 
in e~ N and vN scattering. At low Q^, a reasonable description of the proton and neutron elastic 
form factors is given by the dipole approximation. The dipole approximation is a lowest-order 
attempt to incorporate the non-zero size of the proton into the form factors. The approximation 
assumes that the proton has a simple exponential spatial charge distribution, p{r) = /3oe~''/''°, 
where ro is the scale of the proton radius. Since the form factors are related in the non-relativistic 
limit to the Fourier transform of the charge and magnetic moment distribution, the above p{r) 

yields the dipole form defined by: ^^'^(Q^) = C^.^/J^l + -^-^j . Here C7^'^= (l,^^), OA = 

-1.267, = 0.71 {GeV/cf, and Ma = 1.015 GeV/c^ (see below). 

Since Ma is not equal to My, the distribution of electric and axial charge are different. 
However, the magnetic moment distributions were assumed to have the same spatial dependence 
as the charge distribution {i.e., form factor scaling). Recent measurements from Jefferson Lab 
show that the ratio of ^^^p falls at high challenging the validity of form factor scaling 

and resulting in new updated parameterizations of the form factors |;2j). In this paper we 
present parameterizations that simultaneously satisfy constraints at low where the spatial 
structure of the nucleon is important, and at high where the quark structure is important. A 
violation of form-factor scaling is expected from quark-hadron duality. We use our new vector 
form factors to re-extract updated values of the axial form factor from a re-analysis of previous 
neutrino scattering data on deuterium and present a new parameterization for the axial form 
factor within the framework of quark-hadron duality. 

The new parameterizations presented in this paper are referred to as the duality based 
"BBBA07" parameterization. Our updated parameterizations feature the following: (1) 
Improved functional form that adds an additional dependence using the Nachtman scaling 



variable ^ to relate elastic and inelastic form factors. For elastic scattering (x = 1) ^P'^'^ = 
where Tp,„,iv = Q'^/AMl^^. Here Mp^n,N are the proton (0.9383 GeV/c^), 

neutron (0.9396 GeV/c^), and average nucleon mass (for proton, neutron, and axial form factors, 
respectively). (2) Yield the same values as Arrington and Sick [3J for < 0.64{GeV/c)'^ , while 
satisfying quark-hadron duality constraints at high-Q^. 

For vector form factors our fit functions are Aj\f{(,) (i.e. Aep{(,^), AMpiC'^), AEni^^), 
^Mn(C")) multiplying an updated Kelly[2] type parameterization of one of the proton form 
factors. The Kelly parameterization is: ^^^"^'(Q^) = J2k=o(^kTp /I + J2]^^i bkTj;, where oq = 1 
and m = 1. In our analysis, we use all the datasets used by Kelly[2^, updated to include the 
recent BLAST results, to fit Gep, Geu-, GMp/fJ-p, and GMn/fJ-n {fJ-p = 2.7928, /in = —1.9130). 
Our parameterization employs the published Kelly functional form to G^^^^, and an updated 
set of parameters for G^p'^""^'^((5^). The parameters A]\i{£^) is given by 

n ^ d e 

^Nio = n j^- 

j=l k=l,ky^j ^'^ 

The are equidistant "nodes" on an interval [0, 1] and pj are the fit parameters that have an 
additional property Aj\f{(^j) = pj. The functional form A]\f{S^) (for Gep, Gup, Geu, and Gmu) is 
used with seven pj parameters at ^j=0, 1/6, 1/3, 1/2, 2/3, 5/6, and 1.0. In the fitting procedure 
the parameters of j4jv(^) are constrained to give the same vector form factors as the recent low 
fit of Arrington and Sick |i3j for < 0.6A{GeV/c)'^ (as that analysis includes coulombs 
corrections which modify Gep, and two photon exchange corrections which modify Gmp and 
Gmu)- Our fits to the form factors are: 

GMp{Q')/f^P = AMp{e) ^ gX'"\q^) 
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GEpiQ') = AEpie) X G'^f ^Q') 



GMn{Q')/^^n = Af/:\e)xGMp{Q')/f^p 

GEniQ') = 4r (n X G^p(Q2) X (^^) , 

where we use our updated parameters in the Kelly parameterizations. For Geti the parameters 
a=1.7 and b=3.3 are the same as in the Galster[2^ parametrization and ensure that dGEn/dQ^ 
at for = is in agreement with measurements. The values A{(^)=pi at ^i=0 (Q^ = 0) 
for G^fp, Gep, Geu, Gmu are set to to 1.0. The value A{S)=pi at ^j=l (Q^ — > oo) for Gmp 
and Gep is set to 1.0. The value A{S)=pj at Cj=l for Gmu and Geu are fixed by constraints 
from quark-hadron duality[l].The parameters and plots of the new form factors G_Bp, GAfp//Xp, 
GMn/Mn, and Geu are given in ref.[T]. 

Using our updated BBBA200725 form factors and an updated value qa = -1.267, we perform 
a complete reanalysis of published quasielastic |6, (QE) data on deuterium {v^ n — > //~p) 
using the procedure described in detail in ref. 1^. We extract new values of (given 
in Table [l]), and updated values of FAiQ"^) ■ The average of the corrected measurements of 
Ma from Table |l] is 1.0137 ± 0.0264 GeV/c^. This is to be compared to the average value of 
1.0140 ±0.0160 GeV/c^ extracted from pion electroproduction experiments after corrections for 
hadronic effects, world average value from v^d, V^Ji and pion electroproduction experiments of 
Ma = 1.0137 lb 0.0137 GeV/c^. This is smaller than the recent results[9 from MiniBoone on a 
carbon target {Ma = 1.23 ±0.20 GeV/c^) and K2i^[10: on oxygen {Ma = 1.20 ±0.12 GeV/c^). 
Both collaborations use updated vector form factors, collaborations attribute a difference from 
deuterium to nuclear effects. However, there is experimental and theoretical evidence PJJ that 
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Table 1. M4 [GeVjc?) values published by f^-deuterium experiments [6] and updated 
corrections AM^ when re-extracted with updated BBB 2'o\X\ form factors, and ga=-1.267. 
Also shown is updated Ma from I7^Hydrogen n [7 . 





Figure 1. (a) ^^(Q^) re-extracted from neutrino-deuterium data divided by G^(Q^). (b) 
F^(Q^) from pion electroproduction divided by G^(Q^), corrected for for hadronic effects. Solid 
line - duality based fit; Short-dashed line - Fa{Q'^)a2=V2- The long-dashed line is Fa{Q'^)ai=vi- 
Dashed-dot line - constituent quark model 



Ma in nuclear targets is the same (or smaller) than in deuterium. This discrepancy is important 
for v oscillations experiments since it affects the normalization (at high energies the QE cross 
section is approximately proportional to Ma) and non-linearity of the QE cross section, which 
is relevant to the extraction of v mass difference and mixing angle. 

For deep-inelastic scattering, the vector and axial parts of F2 are equal. Local quark- 
hadron duality at large implies that the axial and vector parts of i^l'"****^ are also equal: 
[Fa{Q^)a2=V2? = {G^)HQ') + rN{Gl{Q^))y{l + r^), where G^Q') = Gep{Q^) - GEn{Q^) 
and G\,{Q^) = GupiQ^) - Gmu^Q^). 

We extract values of Fa{Q'^) from the differential cross sections using the procedure of 



I a) V n ^ u- p 



£ 2 
O 1.8 
1.6 
1.4 
1.2 
1 

0.8 



TTfll 0.05 0.1 



0.5 1 5 



0.6 
0.4 



"T" 



« AVERAGE Deuterium 
-■- M, = 0.8GeV 
- - M, = 1.2GeV 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 2. Average FAiQ"^) /Gd{Ma = 1.015) from (a) i^-deuterium experiments, and (b) 
from pion electroproduction experiments. Also shown are the average E734 and results 
from carbon. Solid line - duality based fit; Short-dashed line - Fa{Q'^)a2=V2'-, long-dashed)- 
Gd{Ma = 1.20)/Gd(Ma = 1.015); dashed-dot line - Gd{Ma = 0.80)/Gz)(Ma = 1.015) 



ref. [1]. The overall normalization is set by the theoretical QE cross section [T2]. We 
then do a duality based fit to i<A(Q^) (including pion electroproduction data) of the form: 

We impose the constraint A^pj^[S^i = 0) = pi = 1.0. We also constrain the fit by 
requiring that A^pj^{^^) yield Fa{Q'^) = Fa{Q^)a2=V2 by including additional "fake" data 
points) for ^ > 0.9 (Q^ > 7.2{GeV/cf). Figure 1(a) shows Fa{Q^) extracted fr om neutrino- 
deuterium experiments divided by G^{Q'^)[lT . Figure 1(b) shows Fa{Q^) extracted from pion 
electroproduction experiments divided by G;^((5^)[12]. These pion electroproduction values 
can be directly compared to the neutrino results because they are multiplied by a factor 
Fa{Q^,Ma = 1.014 GeV/c^)/FA{Q'^,MA = 1.069 GeV/c^) to correct for AMa = 0.055 GeV/c^ 
originating from hadronic effects. Figure (2) shows the same data averaged over all experiments 
in bins of Q^. The solid line is our duality based fit. The short-dashed line is Fa{Q'^)a2=V2- 
The dashed-dot line is a constituent-quark modeljH] prediction. 

In summary, our new parameterizations are useful in modeling v interactions for oscillations 
experiments. Our predictions for GEn{Q'^) and Fa{Q'^) can be tested in future e — N and i^-N 
experiments. 
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